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Overview 


*  Motivation  and  Background 

*  What  is  a  Collisionless  Shock  Wave? 

*  Features  of  the  Collisionless  Shock 

*  The  Shock  Simulation  Model 

*  Simulation  Results  and  Verification 

*  Future  Work 
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Background 


Who  am  I? 

•  Grew  up  in  Olympia,  WA,  The  City  of  Champions 

•  B.S.  Aerospace  Engineering  2016  |  UW,  Seattle 

•  Start  Graduate  Program  at  UW  Autumn  2016 

•  Hobbies: 

-  Studying  classical  guitar, 

-  Observing  the  sky, 

-  Appreciating  existence  and  nature. 
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Investigation  Motivation 


W 


What's  the  Point? 


•  The  ID  Collisionless  Ion-Acoustic  Shock  Wave: 

-  Has  a  simple  premise 

-  Is  rich  in  nonlinear  effects  difficult  for  theoretical  prediction, 

•  Wave  Dispersion 

•  Wave-Particle  Interaction 

•  Various  Wave  Dissipation  Mechanisms 

-  Shock  structure  is  an  active  area  of  research.  In  particular, 

•  Capacity  to  create  ion  beam  of  super-shock  velocity. 

•  Validation  of  Simulation  Frameworks 

-  The  shock  is  a  model  problem  for  simulation  code  validation. 
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Electrostatic  Shock  Features 

*  Propagation  occurs  via 

collective  electric  interactions  lon  Density  Electron  Density 


•  Shock  front  is  dispersed  across 
many  Debye  lengths 

•  Wave  is  accompanied  by  an 
electric  field  perturbation 

•  Two  classes  of  shock: 

-  Subcritical:  M  <  M*  -  1.8 


•  Accompanied  by 
undulating  wave-train 

-  Supercritical:  M  >  M*  -  1.8 

•  Increase  in  dissipation 

•  Reflects  almost  all  ions 


AFRL^I 
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Details  on  Shock  Structure 


Shock  Thickness 


-  Ion-acoustic  dispersion  relation  predicts  the  electron  contribution  to  acoustic  speed 
to  decay  for  wavelengths  approaching  the  Debye  length. 


-  Imagine  shock  formation:  nonlinear  steepening  of  an  acoustic  disturbance. 


Wave  steepening  is  balanced  by  dispersion; 

short  wavelength  components  recede  from 
the  wavefront. 
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Linearized  Ion-Acoustic  Dispersion  Relation 
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©y  Details  on  Shock  Structure 


Wave  Train  Properties 

-  Linearized  Model  of  Dispersive  Shock,  the  Integrated  Airy  Function: 

*  Linearization  of  the  fluid  equations  with  piston-like  boundary  conditions  gives  a 
solution  for  the  shock  behavior. 

•  Assumes  cold  upstream  ions,  therefore  neglecting  shock  reflection. 
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Details  on  Shock  Structure 


Wave  Train  Properties 

-  Linear  Model  Predicts: 

•  The  wave  train  spreads  out  in  time;  peak-to-peak  wavelength  increases. 

•  No  variation  of  non-dimensional  wavelength  with  electron  temperature. 

•  The  wave  train  increases  in  frequency  towards  the  back  of  the  train. 


3-,  Variation  of  First  Ripple  Wavelength  with  Time  in  Linearized  Shock  Profile 
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Details  on  Shock  Physics 


Sources  of  Collisionless  Wave  Dissipation 

-  Landau  Damping: 

*  A  form  of  wave-particle  resonance.  Resonant  particles  sap  wave  energy. 

*  Damping  is  proportional  to  slope  of  distribution  function  at  wave  velocity. 

*  Higher  electron  temperature  means  greater  shock  speed,  so  less  damping. 

-  Ion  Reflection  /  The  Shock  Foot: 

•  Higher  shock  strength  brings  a  greater  potential  difference. 

•  Ions  unable  to  traverse  difference  are  reflected  downstream. 
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Details  on  Shock  Physics 


Phase  Space  Indicators  of  Wave  Dissipation 

(Log  View  of  Distribution) 
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Dissipation  Controls  Wave  Train 


Under-  and  Over-damped  Shocks 


-  Under-damped: 


-  Over-damped: 


Dissipation  is  weak,  ripples  persist. 
High  — 

Ti 


Strong  dissipation  damps  ripples. 


Low 


Te 

Tt 
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The  Simulation  Model 


Initial  Conditions:  Ion-Acoustic  Shock  Tube 

Shown:  Velocity  Distribution  Function  (VDF)  f(x.v) 
And  Potential  Profile  (f>(x) 

High  Particle  Density  Low  Particle  Density 


Kinetic  Model 


Velocity 

A 


Position 


Potential 

A 


Initial  Condition 


Initial  Condition 


Vlasov  Equation 


df  ,  d f  eV <fi  df 


dt 


+  v 


dx 


m 


dv 


=  0 


Poisson  Equation 

Boltzmannian 
Electrons 

Simulation  Procedure: 


e0V2</>  +  p  =  0 


n„  — 


0 

nref  ■  exp  — 


Solve  Poisson's  equation  for  potential. 


-  VDF  for  ion  density 


<t> 


-  Electron  density  from  exp  — 

-L  e 

Update  the  distribution  function  via 
Vlasov  equation. 


Position  • 
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Repeat. 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 

Density  Difference:  30% 

AF  RL  13 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 

Density  Difference:  30% 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 

Density  Difference:  30% 

AF  RL  15 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 

Density  Difference:  30% 

AF  RL  16 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 

Density  Difference:  30% 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 

Density  Difference:  30% 

AF  RL  18 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 

Density  Difference:  30% 

AF  RL  19 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 
Density  Difference:  30% 


. . . 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 

Density  Difference:  30% 
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Typical  Simulation  Result 


Lightly  Damped  Case 

Electron  Temperature:  5  eV  |  (Te/Ti  =  25) 

Density  Difference:  30% 
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Typical  Simulation  Result 
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Heavily  Damped  Case 

Electron  Temperature:  1.5  [eV]  |  (Te/Ti  =7.5) 

Density  Difference:  30% 
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Typical  Simulation  Result 


1.32e+15-i 

1.3e+15- 

1.28e+15- 

1 .26e+l  5  - 

1 .24e+15- 

1 .22e+l  5  - 

1.2e+15 

1.18e+15 

1.16e+15 

1.14e+15 

1.12e+15 

l.le+15 

1.08e+15- 

1.06e+15- 

1.04e+15- 

I.02e+I5- 

le+15- 

v.8e+i4  — 

-0.12 


A 


L 


Heavily  Damped  Case 

Electron  Temperature:  1.5  [eV]  |  (Te/Ti  =7.5) 
Density  Difference:  30% 
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Typical  Simulation  Result 


Heavily  Damped  Case 

Electron  Temperature:  1.5  [eV]  |  (Te/Ti  =7.5) 

Density  Difference:  30% 
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Model  Verification 


Comparison  with  Linearized  Solution 


-  Evolution  of  the  First  Ripple  Wavelength: 

•  Simulated  weak  shocks  (1%  jump)  for  different  temperature  ratios. 

•  Greater  ion  reflection  at  lower  electron  temperature  alters  the  first  ripple. 

•  Agreement  improves  with  hot  electrons,  as  ion  reflection  decreases. 
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Model  Verification 


Comparison  with  Linearized  Solution 


-  Evolution  of  the  First  Ripple  Wavelength: 

•  Stronger  shocks  (4%)  show  poorer  agreement  at  lower  electron  temperature. 

•  Stronger  shocks  reflect  more  ions,  causing  stronger  dissipation. 

•  Agreement  with  linearized  model  improves  with  electron  temperature. 
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Model  Verification 


Comparison  with  Linearized  Solution 

-  Trailing  Ripple  Wavelengths: 

•  Subsequent  wave  train  ripples  diverge  from  the  model. 

*  Simulated  ripples  do  not  approach  zero  wavelength  like  the  linearized  model. 
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USCS 1  |  Mach  Number 


Results  on  Shock  Structure 
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Other  Weak  Shock  Results 


-  Shock  Mach  Number: 

•  For  the  same  density  jump  (1%  here), 

-  Shock  speed  increases  with  electron  temperature,  Mach  number  decreases. 

•  The  first  wave  train  oscillation  peak  travels  less  quickly  than  the  shock  ramp. 
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USCS  1  |  Mach  Number 


Results  on  Shock  Structure 


Other  Weak  Shock  Results 


-  Shock  Mach  Number: 

•  The  same  trends  are  apparent  for  slightly  stronger  shocks  (4%  jump  here). 

*  Ion  reflection  appears  to  damp  the  shock  acceleration. 
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Model  Validation 


Gathering  Experiment  Data 


-  Double  Plasma  Device: 

*  Can  produce  low  Mach  number  collisionless  shocks. 

-  Diagnostics: 

*  Langmuir  probes;  density  and  electron  temperature. 

-  Shock  speed,  and  density  at  fixed  position. 

*  Ion  temperature  is  necessary  for  proper  validation. 

Validation  Comparisons 


-  The  Wave  Train: 


*  Simulation  suggests  ripples  should  occur  for: 

-  Moderate  Mach  numbers  (~1.2  -  -1.6) 

-  High  electron-ion  temperature  ratio  (>10) 
-  Wave  Train  Wavelength 


Density 


Simulation 


-  Shock-Front  Mach  Number 

-  Reflected  Ion  Beam  Velocity 


Time 


Experiment 
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5  '  Towards  Validation  and  Extension 


Towards  Validation  of  Collisionless  Model 

-  Verification: 

•  Simulated  weak  shocks  agree  with  the  linearized  cold-ion  solution. 

-  Validation: 

•  Results  should  be  compared  to  experimental  shock  data. 

•  Inconsistencies  in  published  1969  double-plasma  device  data  hampered  validation. 

Future  Work:  Extension  to  Moderately  Collisional  Problem 

-  Collision  Operators: 

•  Presence  of  a  few  mean  collisions  in  shock  transition  invalidates  Vlasov  equation. 

•  The  collisional  ion-acoustic  shock  requires  the  full  Boltzmann  equation  with 
appropriate  choice  of  collision  operator. 

•  Moderately  collisional  shocks  should  be  hybrid  dispersive-dissipative  shocks. 

-  Validation  of  Collision  Algorithms: 

•  Tuning  parameters  in  a  shock  experiment  can  yield  moderately  collisional  data. 
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Results  on  Shock  Structure 


Example  of  Shock  Speed  at  1% 


I 
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Details  on  Shock  Structure 


Wave  Train  Properties 

-  Sagdeev  Pseudopotential: 

*  An  analogy  to  motion  in  a  potential  well,  describes  oscillations  in  wave  train. 

*  Dissipation  within  the  well  leads  potential  to  settle  at  the  minimum. 
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Electric  Field  Details 


fc  |  Wavenumber  [^.]  le3 
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